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We investigate the stability of D dimensional singly rotating Myers-Perry-AdS black holes under 
superradiance against scalar field perturbations. It is well known that small four dimensional ro¬ 
tating or charged Anti-de Sitter (AdS) black holes are unstable against superradiance instability of 
a scalar field. Recent works extended the existence of this instability to five dimensional rotating 
charged AdS black holes or static charged AdS black holes in arbitrary dimensions. In this work we 
analytically prove that, rotating small AdS black holes in arbitrary dimensions also show superra¬ 
diance instability irrespective of the value of the (positive) angular momentum quantum number. 
To do this we solve the Klein-Gordon equation in the slow rotation, low frequency limit. By using 
the asymptotic matching technique, we are able to calculate the real and imaginary parts of the 
correction terms to the frequency of the scalar field due to the presence of the black hole, confirming 
the presence of superradiance instability. We see that, unlike in the case of static AdS black holes, 
the analytical method is valid for rotating AdS black holes for any value of angular momentum 
number and space-time dimensions. For comparison we derive the corresponding correction terms 
for Myers-Perry black holes in the black hole bomb formalism in Appendix and see that the results 
are in agreement. 

PACS numbers: 04.20.Jb; 04.40.Nr; 11.27.+d, 04.50.-h 


ozgur.delice@marmara.edu.tr 



2 


I. INTRODUCTION 


Superradiance is a phenomenon where a field satisfying certain conditions is amplified by its interaction with a 
dissipative system [1], In the black hole superradiance [2, 3], the interaction of scalar, electromagnetic or gravitational 
fields with the black hole event horizon, which behaves as a one way membrane mimicking a dissipative system, 
causes the field to be amplified, if the field satisfies the superradiance condition. The condition for superradiance for 
Kerr black holes is u> < mflh [4-6] where ui is the frequency of the wave, m is the azimuthal number, and fi/, is the 
angular velocity of the horizon of the black hole. A possible indirect observation of superradiant scattering from black 
hole-pulsar binary systems were proposed recently [7]. The physics of superradiance and several aspects of black hole 
superradiance is discussed in an excellent review [1] and we refer this review for further information. 

The superradiance instability requires a mechanism such as a potential barrier, which localizes the field near the 
horizon of the black hole and does not permit scattered waves to escape to infinity. The continuous reflection of 
the amplified waves from the barrier towards horizon may lead to an instability of the black hole, i.e., its rotational 
energy and angular momentum is decreased by such a process. It is well known that the Schwarzscliild black holes are 
stable [8, 9] against such perturbations. More general black holes, however, may develop instabilities under certain 
conditions. For example, although the Kerr black hole is stable under massless scalar, electromagnetic or gravitational 
perturbations [10, 11], a massive scalar field may cause superradiant instability on this black hole [12-16], since the 
mass term acts as a natural mirror. In order to imitate the conditions leading to the superradiance instability, the 
black hole bomb mechanism is often used [17-33]. In that mechanism the black hole is though to be surrounded by a 
hypothetical reflective mirror leading the field to be amplified back by the black hole, resulting an instability. The use 
of this hypothetical setting helps to analyze certain aspects of this instability such as the time scale of the instability 
and the relation of the instability with the parameters of the black hole, the field or the mirror. 

Although their astrophysical importance is limited, the black holes living in an AdS spacetime have become very 
important after the gauge/gravity duality is introduced [34]. Thus, it has become a necessity to explore the stability 
of these black holes under different perturbations, such as the superradiance instability. Since the boundary of AdS 
spacetime behaves as a mirror wall, a black hole living in an AdS spacetime may be unstable against superradiance 
as well, similar to the black hole bomb mechanism. Actually, it is well known that large AdS black holes are stable 
[35]. However, the four dimensional rotating small [36, 37] or charged [38] AdS black holes are unstable against 
superradiance instability of a scalar or charged scalar field. In five dimensions, it was shown that small charged AdS 
black holes with two rotation parameters are also unstable [39]. This instability is also reported for five dimensional 
hairy AdS black holes in [40]. Superradiance also leads to gravitational instabilities for AdS black holes in four and 
higher dimensions [41-46]. Recently, small charged static AdS black holes in D dimensions are also shown to be 
unstable [47] against superradiance, together with the correction of an erroneous conclusion of [39] that for certain 
values of orbital number l, the instability is not triggered in five dimensions. Actually, this is corrected using a 
numerical approach, since for those particular values of l, analytical methods fail for static black holes as in their 
case. See also [28] for five dimensional analytical treatment of such values of l in the black hole bomb mechanism. In 
this paper, our aim is to investigate the stability of rotating AdS black holes in arbitrary dimensions against scalar 
perturbations using the analytical methods. We will especially show that, for rotating black holes, unlike static ones, 
the analytical method is applicable for any (positive) value of orbital quantum number. We also compare these results 
with the superradiance instability of D dimensional singly rotating Myers-Perry black holes under scalar perturbations 
in the black hole bomb mechanism which we review in the Appendix and find a perfect agreement. 


II. SUPERRADIANCE INSTABILITY OF MYERS-PERRY-ADS BLACK HOLES 

A. Metric and its properties 


A general D = 4 + n dimensional rotating-AdS black hole is given by the Myers-Perry-AdS solution [48-50]. Here 
we consider the special case where only a single nonvanishing rotation parameter a [48, 51] exists. This spacetime is 
described by 
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where M is the mass parameter, £ = \J—(D — 2 ){D — 1)/(2A) is the AdS radius, A is a negative cosmological constant, 
and dfl 2 is the standard metric of the n dimensional unit sphere. The metric functions are given by 


E = r 2 + a 2 cos 2 0, A r = (r 2 + a 2 ) (\ + ^ - 2Mr 1 ””, 
a 2 2 _ a 2 

Ae = 1 ~-t 2 C0S 


( 2 ) 

( 3 ) 


This black hole has physical mass /z and angular momentum J given by [52, 53] 
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where 4^-2 is the surface area of a unit D — 2 sphere. The angular velocity of the horizon with respect to rotating 
infinity is given by 


fi/t = 


( 5 ) 


and corresponding expression for the nonrotating infinity relevant for black hole thermodynamics is 12 = + a/£ 2 . 

The event horizon of this black hole, r = r> i; is located at the largest root of the equation A r (r) = 0. There is a 
Bogomornyi-Prasad-Sommerfield (BPS)-like upper bound [54] on the rotation parameter, i.e. |a| < £, otherwise the 
metric describes a naked singularity. These black holes also suffer from an ultraspinning instability [55, 56]. It was 
demonstrated in [41] for four dimensional Kerr-AdS black holes that superradiance instability is present if the horizon 
angular velocity satisfies fl£ < 1 and the end point of the instability is described by a Kerr-AdS black hole whose 
boundary is an Einstein universe rotating with the speed of light or equivalently the corotating Killing vector becomes 
space-like and in equilibrium with a scalar field cloud. 


B. Klein-Gordon equation 

In this subsection we consider the Klein-Gordon equation for a scalar field 4> with mass m of the form 

-m 2 4> = 0. (6) 

Since this equation is known to be separable [48, 57, 58] for general Myers-Perry-AdS black holes, we can consider 
the following ansatz for the scalar field 


$ = e irn< t>~ iult Y (fi)0(0)i?(r). 

The resulting equation separates into its angular and radial parts as follows: 
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Here the terms j (j + n — 1) are the eigenvalues of the spheroidal equation for Y of the n-sphere [59] with j assuming 
integer values. We need the eigenvalues \ji m of the angular equation (8) where the eigenvalues of this equation for 
four [60], five [39] and higher dimensions [61] is discussed recently. In the slow rotation low frequency limit, i. e. 
auj <C 1, and for the case the mass of the scalar field is very small such that m 2 a 2 ~ 0, the eigenvalues Xji m of the 
Kerr-AdS spheroidal harmonics can be expanded into a Taylor series as 


OO 
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where u> = s/to 2 — m 2 , l is a positive integer which satisfies l > j + |m|, and certain terms of f p are given explicitly 
in [60, 61]. We take m as a positive integer in this paper. The crucial point here is that Xji m have small correction 
terms added to l terms, which prevent the eigenvalues from being an exact integer, unless the rotation parameter a 
vanishes, as in the case of static black holes. The radial equation cannot be solved analytically but using the matched 
asymptotic expansion method, presented first by Starobinsky [5], we can obtain the asymptotic solutions near the 
black hole horizon and in the far region. Their matching at the intermediate region will enable us to calculate the 
correction term to the frequency of the wave and subsequently determine the existence of the instability. 


C. Asymptotic behaviour of the scalar held and the superradiance condition 

In order to discuss the asymptotic behavior of the scalar field, it is useful to introduce a tortoise coordinate r* by 
considering the transformation r = r(r*) of the radial coordinate and the wave function as follows: 

7 Z(r) = \Jr n (r 2 + a 2 ) R, (11) 

/7r** r*^ _i_ q 2 

= (12> 

These transformations bring the radial equation into a Sclrrodinger-like form, 

^ + V(r)H = 0, (13) 


where the potential term V (r) is given by 
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Near the horizon r —> rh, since A r ~ 0, the potential term becomes 

V(r-> r h ) ~ (w - mVL h ) 2 . (15) 

This implies that the solution near horizon becomes, in the classical limit where only the ingoing waves present, 
R(r —> rh) ~ e -iut-z(u-mSi h )r anc j ^ ie superradiance condition is 

u> < milk, (16) 

where the field is amplified when the frequency of the wave satisfies this condition. 

As r —> oo the potential becomes infinitely large, implying the vanishing of the scalar field at radial infinity, i.e. 

$(r -A oo) -a 0. (17) 

As a consequence of these observations, the appropriate boundary conditions turn out to be the Dirichlet boundary 
condition at the radial infinity and incoming wave boundary conditions on the horizon of the black hole. 


D. Near region solution 

Here we solve the radial part of the Klein-Gordon equation at the region near the event horizon of the black hole 
r ~ rh in the slow rotation low frequency limit, noting that near the horizon the effects of the cosmological constant 
are negligible. We also assume that the Compton wavelength of the perturbations are much larger than the size of 
the horizon, i. e., mr> ( -C 1 . In these limits, i. e., r-ri,«l,wa~0,a 2 ~0,ri,<<,«« l, Eq. (9) becomes 

r-"! ^ (w - mn h f - C R(r) = 0, (18) 

where 

C = Xji m + 1 T^ r h ~ l{l + n + 1) + e. (19) 
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Note that the term e carries all the correction terms due to the rotation of the black hole and mass of the scalar field. 
This small term is not arbitrary and can be calculated using (10) for given values of the parameters of the black hole 
and the scalar field. We will keep this term since it will be crucial in the forthcoming analysis. In order to obtain 
solutions of Eq. (18), let us consider a new radial coordinate 

x = r n+1 , (20) 


which brings the radial equation into 


(n + 1) 2 A-^- 

ax 



2(rt + 2) 

a: + " +1 (cj — mflh) 2 



R(r) = 0 , 


( 21 ) 


where A = r 2n A r w x 2 — 2Mx + a 2 x 2n ^ 1 + n ) = (x — x+)(x — x-). In four dimensions n = 0, we have two horizons, and 
in the higher dimensions we can set the last term in the approximation and the inner horizon to zero. Nevertheless, 
in order to have a unified treatment we keep the X- term in our expressions for higher dimensions as well. To bring 
Eq. (21) into the form of a hypergeometric equation, we first define a new dimensionless variable 


x — x + 
x — X- 


( 22 ) 


with the event horizon being at z = 0. Then the radial equation becomes 
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where 
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Finally, if we define the radial part of the wave function R in terms of a new function R(z) as follows 


1 + V1+4C 

R = z lu (l-z)—i — R(z), 


(26) 


we see that the radial part of the Klein-Gordon equation (23) is equivalent to a hypergeometric differential equation 
given below: 

z(! - z )~j^ + [7 - (« + /? + !) z] ^ - a/3T = 0, (27) 

where the constant parameters a, and 7 are given by 


a = 
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1 + sjl + 4C 
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+ 2 i cj, 


7 = 1 + 2iCj. 


(28) 


The most general solution of this equation in the neighborhood of z = 0 is given by 

T(z) = Az 1 “ 7 F(a -7 + l,^- 7 +l, 2 - 7 ,z) + BF(a,ft 7 ,z), (29) 

where A and B are arbitrary integration constants and F(a, j3,j, z) = 2 -fi(a, /3; 7 ; z) is an ordinary hypergeometric 
function [62]. By considering (26), the general solution of the radial Klein-Gordon equation near the horizon can be 
written as 


1 + V1+4C 

R = Az lUJ (l — z) 2 F(a — 7 + 1,/? — 7 + 1,2 — 7 , z) 
+ B z iQ (l- z) 11 ^ F(a,p n ,z). 
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The first and second terms represent an ingoing and outgoing wave respectively at the horizon z = 0. Because we are 
working at the classical level, there cannot be outgoing flux at the horizon. Therefore the second term of the solution 
should vanish, which requires B = 0 and the radial part of the solution near horizon becomes 


R = Az~ ia, ( 1 - z) 1 +V 2 1+4C p 


1 + V1 + 4C 1 + V1 + 4C 


— 2iui, 1 — 2 zw, z 


(31) 


Since we will match the near and the far region solutions at the intermediate region, we need the large r, i. e., 
z —> 1 limit of this solution. We can use the z —>■ 1 — z transformation law [62] for the hypergeometric functions 


F(a, 6 ,c, z) = ° —7T f («: b,a + b- c+1,1- z) 

I (c — a) 1 (c — b) 

+ (i - g)e _ a _ b r(c)r(a + 6-c) f(c _ ac _ bc _ a _ b + ll _ z) 


(32) 


Note that this formula does not work for c — a — b = ±fc, where k is an integer and for this case a transformation 
involving logarithmic terms is required [62], which will make the matching with the far solution at intermediate regions 

impossible. For our case, however, since a = (1 + \J 1 + 4£)/2, b = (1 + \J 1 + 4£)/2 — 2 iCj and c = 1 — 2 iui, we have 

c — a — b = — \Jl + 4£. Since for any value of l, £ is not an integer, contrary to the previous claim [47], unlike static 
black holes, the above transformation holds for any value of l for rotating black holes. Thus, we do not end up with 
any logarithmic terms in the above transformation. We can also see this if we replace £ = 1(1 + n + 1) + e in the above 
expression; we find that 


c — 


a —b = —y 1 + 4£ 



2(1 + e) 
n + 1 


+ 0 (e 2 ), 


(33) 


where 


e 

2 / + n + 1 


(34) 


The resulting expression can be an integer if and only if e vanishes and 1/(1 + n) is half integer. Here e is the term 
denoting small corrections to the eigenvalues of angular equation due to the rotation of the black hole and mass of 
the scalar field given in (10) and (19). It is not an arbitrarily small parameter, its value can be calculated in principle 
using (10). In summary, unlike static black holes, we can use the transformation law (32) for any value of l, keeping in 
mind that there are always small correction terms to be added to the integer values l. making the relevant expressions 
close to an integer up to a small correction term. 

Thus using the above transformation, by using the property that F(a, b, c, 0) = 1, and considering that the Compton 
wavelength of the scalar field should be much larger than the black hole horizon implying mr b <C 1 , the large r limit 
of the near region solution can be expressed as follows: 


with 


R(r) = A i r l +A 2 r~ l - n - 1 , 
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ip - p 1 

where we have also used the fact that 1 — z —> + r i + n ~ — for r —> oo to derive the above expressions. For clarity we 
discarded e terms in the above expressions except for the terms which might have poles if e vanishing. Actually, as 
we will see, the erroneous conclusion in [39] is not caused by the poles in the above equations but not taking these 
small corrections into the last step of the calculations. 
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E. Far region solution 

For the far region r-fk >r/,,we can ignore the effects of the black hole by setting a = M = 0, then the radial 
equation reduces to [47]: 


r 2 \ d 2 R 
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r n + 2] dR 
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geometric differential equation of the form (27) with 
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The regular solution of this equation, obeying the boundary condition at r = oo is given by 

R(y) = C{ 1 - y) l / 2 y“ e ' 2 ~ a F[a, a - 7 + 1, a - 0 + 1, iT 1 ], (42) 

where C is an integration constant. The small r limit of this solution can be obtained by - —> 1 — y transformation 
of hypergeometric functions [62], which yields 

R(r) = C\r l + C 2 r- l ~ n - 1 , (43) 


c(-i)‘/ 2 r i + sjm 2 £ 2 + (=^) 2 r[-/- 2 ±i] 

^1 = f / I f / - | > ( 44 ) 
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Note that this solution (43) is for pure AdS spacetime for r € [0,oo). The regularity of this solution at the origin 
requires the term to vanish, which can be possible only if one of the T functions at the denominators of C 2 

is equal to oo. This yields a discrete (positive) spectrum [47] 

cun£ = 2 N + n 2 3 + 1 + \jrh 2 £ 2 + 3 ) ’ (46) 

where N = 0,1, 2,.. . is the radial overtone number. In the presence of the black hole horizon, it is natural to expect 
that the frequency of the waves has small (complex) modifications, compatible with the limits we have considered in 
the present analysis, as follows: 

u> = lun + id. (47) 

Replacing this into C 2 and considering the fact that T (—N — i5£/2)~ l = (—l ) w+1 iS£N\/ 2, we have 

C(—1 / /2 (—l) N+1 iSN\ T 1 + ^Jm 2 £ 2 + (^j 2 T[/ + 2±i] 

C 2 = -r-^- J 1 -■ (48) 

2l-(l+n+2) T 2 H +i + jV + y'ro 2 | 2 + (”13 ) 2 



F. Instability 


The near and far region solutions (35) and (43) can be matched at the intermediate region <C r — ?+ <C 1/w, since 
the asymptotic forms of these solutions have the same powers of r. Considering the condition C\/A\ = C 2 /A 2 = 1, 
the matching yields the following expression for the correction term of the frequency 
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This main result is, apart from notational differences and the fact that the superradiance factors are different, in 
agreement with the corresponding expression for the charged static AdS black holes in arbitrary dimensions [47]. 
Although one can make a numerical analysis to investigate the presence of instability using the above expression, we 
can further expand the Gamma functions to prove it analytically. It turns out that the last multiplicative factor in 
the above expression is crucial to determine the real and imaginary parts of the correction to the frequency. Actually, 
we can use the expansion of the expression first presented by us in [27] as 
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r (-uh-2 iu) 
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(50) 


Note that the imaginary part of above equation vanishes if e = 0 and 1/(1 + n) is half integer. This is the reason 
leading to the erroneous conclusion in [39] that in five dimensions, i.e. n = 1, inodes corresponding to odd values of l 
do not trigger the instability. However, since e is not vanishing, the imaginary part of the above equation is also not 
vanishing, which corrects this erroneous conclusion. 

We can evaluate the other gamma functions as given in [47]: 
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Using these formulas, we can express the correction term as a sum of its real and imaginary parts as follows: 

'(/ + e)7rl 


5 = —a -|sinh (27rcD) — i tan 
where the positive multiplicative factor cr has the expression 
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Thus the total frequency of the scalar field becomes 

u> = ujn + iS = Wjv — cr tan 


(l + e)i 
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where ujn is given in (46). The real and imaginary parts of the frequency hence become 

'{l + e)n] 


94c[w] = ojn — <J tan 


1 + n 


cosh (27rw), 


3m[w] = —asinh (27rd)). 


(56) 

(57) 


The results presented above confirm that, irrespective of the value of l , the rotating AdS black holes always suffer 
from the superradiance instability whenever the frequency of the scalar field satisfies the superradiance condition 
ui < 0 which means u> — m flh < 0, since we have 


<j> ^ e~ iuJt = e 3m M t 


(58) 


Thus, the imaginary part of the frequency is always positive under this condition, resulting in an exponential growth 
of the field. The time scale of this exponential growth is inversely proportional to the imaginary part of the frequency, 
i. e., t ~ l/3m[u;]. This fact allows us to obtain a relation between the AdS radius of the spacetime and the instability 
time scale as 


r ~ e 2 (J+1)+n . (59) 

Hence 3m[w] decreases and r increases with increasing n for fixed AdS radius. This implies that as D —> oo, r —> oo 
and the instability becomes ineffective in the large D limit. 

The onset of the instability is determined by 9fe[w] = rnilh. Since the real part of the frequency is inversely propor¬ 
tional to the AdS radius, there is a critical value of t'o such that a superradiant wave starts becoming nonsuperradiant. 
This particular value is 

If we compare our results with the black hole bomb mechanism of Myers-Perry black holes [19, 26], which is rederived 
in the Appendix for a comparison, we see that there is a perfect agreement with the corresponding expressions of real 
and imaginary parts of the frequency correction terms |[(56) and (57)1 and [(111) and (112)1) and the time scale of 
the instability [(59) and (113)]. 

Note that, although the above results will not change, some of the gamma functions in the above expressions can 
be further expanded when the term 1/(1 + n) is integer or half integer. Since we have taken care of the possible poles 
in the gamma functions, we do not have to consider small correction terms to the l below. 


1. The case j4 — = k where k is an integer 


For this case, it is possible to further expand the general equation since some of the gamma functions can be easily 
expanded and also we have the following relation: 
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The resulting expression for the frequency is 
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) 2_ 


k N / l \ 

xn^+^viH'+^+i)- 

p=i j =i v ^ 


(60) 


(61) 


(62) 


It is clear that the imaginary part of the frequency becomes positive under superradiance, implying the superradiance 
instability for small rotating AdS black holes. These results and the correction terms are in accordance with the 
previous works [25, 28, 36, 47]. The real correction term is very small compared to ujn and thus can be ignored. 
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2 . 


The case y2_: = k + 1 where k 


is an integer 


For this case we have 


1 + n 


A 1 — 2'iw 


COSh(27Tw] 


k+1 

n 

p =1 


P - 2 ) + 4 ^ 2 


and also 


'( 1 + lTn) = 2 _fc_ 1 v / ^( 2fc + 1 ) !! - 


Using these expressions and expanding the trivial gamma functions, we have 

(1 + n)a' 


UJ = U) N + 


— ia' tanh (27 tcD) 


with 


a = 


7 T (r| +n - r 1 + n ) 2+2k [(2k +1)!!] 2 T 

n±3.+l + N+ \J rh 2 £ 2 + (2±^) 2 


2 2 (i+fc) pi+n+2 Nl ( 2 fc + l)!(2fc + 2)! r ^ + 2±1 j r 

1+N + ^m 2 £ 2 + (^) 2 


fe+i 

><n 

p =i 


P - 2 1 + 


TV 


xIK'+^+ji- 

i=i 


(63) 


(64) 


(65) 


( 66 ) 


This result shows that when —^ is an half integer, which can only be possible when spacetime dimensions are odd, 
there is a superradiance instability when the wave satisfies the superradiance condition, because the imaginary part 
of the frequency becomes positive. This result analytically corrects the erroneous claim in [39] that for rotating 
AdS black holes in five dimensions there is no instability for even values of l. Hence, the above expressions for 
D = 5 (n = 1) show that there indeed exists an imaginary part of the frequency which is erroneously claimed to be 
vanishing in [39]. This claim was actually corrected in [47] for charged static AdS black holes using only a numerical 
method, since the analytical method fails for those values of l for static black holes. Thus, our paper completes the 
picture by showing analytically that the instability exists for rotating AdS black holes with the help of the fact that 
the analytical methods are valid for any value of l for rotating black holes. 

Since the real correction term, the second term in (65), involves a divergent term (1 + n) for large n, we need to 
compare its magnitude with lon given in (46) and also its behavior for large values of n. To analyze the ratio of the 
real correction term to normal modes, we consider the case when m = k = r_ = 0 and N = 1 for simplicity, then we 
have 


(1 + n)a l /r + \ 2 ( ra+1 ) (1 + n) 2 (2 + n)T[4 + l + n] 
eiruiN VT/ 2 e (11 + 3n) T[1 + n] T 


(67) 


Note that since the small parameter e is not arbitrary and (r+/£) 2 ( n+1 ' > /e is finite, this ratio is finite and in the 
limits we consider where r+ <C £, the term (r+ /£) 2 ( rl + 1 ) decays exponentially for large n, making the combination 
small. Hence the correction term is several orders of magnitude smaller than normal modes ui n, for example in five 
dimensions (n = 1), for reasonable values £ = 1, r+ = 0.01 and a = 0.001, we have wjv = 7 and the correction term 
is 0.038, which is 0.5% of u>n, where we have chosen e = (owat) 2 . Our analysis also shows that when k or n increases 
this term becomes much smaller compared to wjv, keeping in mind that u> n£ = 21V + n + 3 + l for massless case (46). 
Also for large values of n, the above term gets much smaller and for the n —> oo limit it converges to zero. 


III. CONCLUSION 

In this paper we have investigated the superradiance instability of a scalar field for singly rotating small AdS 
black holes in arbitrary dimensions analytically. This instability is originated from the existence of the superradiance 
mechanism of black holes together with the reflective boundary conditions of the asymptotic infinity of the AdS 
spacetime. As a consequence of this behavior, the boundary of the AdS spacetime behaves like an infinite potential 
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barrier, localizing the scalar field near the horizon leading to an instability. The end point of this instability is expected 
such that we have a rotating hairy black hole having less energy and angular momentum, living in an AdS spacetime 
whose boundary rotating with the speed of light. Thus, the next direction on this research topic might be to explore 
the existence and the properties of these scalar hairy AdS black holes. 

We have analytically shown that for rotating small AdS black holes perturbed by a scalar field satisfying the 
superradiance condition, there is always superradiant instability, irrespective of the value of spacetime dimensions 
D > 4 and the value of orbital number l > 0. Our result generalizes the superradiance instability of a scalar field 
for rotating AdS black holes to arbitrary dimensions. Actually, in the pioneering work [36] proving the instability 
of four dimensional small rotating AdS black holes under scalar perturbations, it was claimed that such instability 
must be present for higher dimensions as well. However, an analytical analysis for rotating small AdS black holes 
has not been presented, except for D = 5. For charged static AdS black holes, the instability is proved for generic 
dimensions [47] together with the observation that, unlike for the rotating AdS black holes as we have shown in our 
paper here, the analytical method fails for certain values of the orbital quantum number. The time scale of the 
instability is proportional to a power of the radius of the AdS spacetime increasing with spacetime dimensions, similar 
to the case of the black hole bomb mechanism. These results are in accordance with the instability of rotating black 
holes surrounded by a hypothetical reflective mirror in the black hole bomb mechanism, in which we have reviewed at 
the Appendix. We have also recovered in the Appendix the previous results that the bosonic and fermionic thermal 
factors for a scalar field absorption naturally arise for rotating black holes for D = 5. 
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APPENDIX: BLACK HOLE BOMB MECHANISM FOR MYERS-PERRY BLACK HOLES 


A. Myers-Perry Black Hole Spacetime 


In order to compare the results we have found for rotating AdS black holes in D dimensions with the rotating 
black holes in black hole bomb formalism, here we review the black hole bomb mechanism for D dimensional rotating 
Myers-Perry black holes, which were studied before in four [19] and generic dimensions [26] and shown to suffer 
from superradiance instability. Since the existence of the notational differences and the real and imaginary parts of 
the frequency correction terms are not calculated explicitly in [26], here we have preferred to derive all the relevant 
expressions in detail. Hence here there are no new results except the explicit expressions of frequency correction 
terms, and the relations of decay rates with bosonic and fermionic thermal factors. This part is partially based on 
[27] with some corrections. We consider the D = 4 + n dimensional Myers-Perry black hole with mass M and a single 
rotation parameter a given by the metric 


2 A — a 2 sin 2 6 2 2a(r 2 + a 2 — A) . , (r 2 + a 2 ) 2 - A a 2 sin 2 6 . 2 2 

as =-—- at -—-sin ddtdcp H-—-sin 6 d<p 


+ —dr z + T,d6 z + r z cos z 0dSl z 
A 


( 68 ) 


where the metric functions are given by 

£ = r 2 + a 2 cos 2 9 , 

A = r 2 +a 2 — 2Mr l ~ n , (69) 


and here again dfl 2 is the standard metric of the ?i-sphere. This black hole has physical mass /i and angular momentum 
J which can be found by setting H = 1 in (4) for the same parameters for AdS black holes. The event horizon of the 
black hole is located on the largest real root of A = 0 given by 


r 2 h + a 2 - 2Mr] l ~ n = 0. 


(70) 


In four and five dimensions the existence of the horizon sets an upper bound (extremal limit) on the rotation 
parameter a (a < M for D = 4 and and a < y/2M for D = 5). In these limits in four dimensions there are two 
horizons r± where they are called the outer and the inner horizons. In dimensions greater than 4, however, there is 
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always one horizon (r = rh) and they can have an arbitrarily high rotation parameter a if D > 5. They are called 
ultraspinning black holes and they have a regular horizon and arbitrarily large angular momentum per unit mass. 
These ultraspinning black holes are known to suffer from several instabilities [63, 64]. In this paper we consider the 
slow rotation regime, in which the ultraspinning instability is not present. 

The temperature of this black hole is given by [65] 


n = r + ~ t (—) , 

4tt V + a r h J 47 T \ r h J 

where the last expression is valid in the slow rotation limit aCl. 


(71) 


B. Klein-Gordon equation 

Here we calculate the massless (m = 0) Klein-Gordon equation for a scalar field $ given by 

= 0 . 


(72) 


If we use the standard separation ansatz 


$ = e im ^- iajt F(fl)e(6»)i?(r), 


we see that the Klein-Gordon equation is separated into angular and radial equations given by 


1 


cos 11 9 sin 0 dO 


cos" 9 sin ) + 
ad 


2 2 2 /i 

a uj cos 6 — 


sin 


j(j + n- 1) 
cos 2 9 


+ A 


jlm 


0 = 0 , 


(73) 


(74) 


1_d_ 

r n dr 



(a 2 + r 2 ) 2 

A 



ma 

a 2 + r 2 


2 


j(j + n- 1 )a 2 



R(r) = 0, 


where the constants Xjim are given by 


Xjim — Xjiui H - a uj 2amcj. 


(75) 


(76) 


As in the AdS case, the terms j(j + n — 1) are eigenvalues of the spheroidal equation for Y of the n-sphere [59] where 
j is an integer. For this case, the angular equation (74) takes the form discussed in [66] and in the slow rotation limit 
we can expand the eigenvalues in a Taylor series of the following form 


A jim — 1(1 + n + 1) + f p (auj) p , (77) 

p= i 


and its first several values are explicitly calculated in [66]. 

We can also express the radial equation in a Schrodinger-like form, 

^ + V(r)K = 0, (78) 

where transformation equations and the potential term are given in (11), (12), and (14) by setting S = 1, rh = 0 and 
replacing A r with A given in (69). Near the horizon, the behavior of the scalar field is the same with the AdS case 
(15), i. e., V(r —> rh) ~ (u> — in^lh) 2 with asymptotic behavior of the field R(r —> rh) ~ e -tut-i(u>-mCi h )r ^ w j iere the 
angular velocity of the horizon for this case is given by 


f \ 


-.2 


(79) 


The asymptotic form of the radial function at radial infinity is different from the AdS case, since, for r —> oo the 
potential becomes V(r —> oo) = w 2 , which yields the solution to be of the form R ~ . it is easy to verify that 

when the well-known condition, 


ui — rnYLh < 0 


(80) 
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is satisfied the waves reflected back from the black hole with increased energy and the phenomena called superradiance 
occurs. The question is whether it is possible to localize waves satisfying this condition near the horizon to lead to an 
instability on the black hole. One of the methods for localizing these waves is considering an artificial reflecting mirror 
surrounding the black hole, i. e., black hole bomb mechanism, in which the waves are continuously scattered between 
the black hole horizon and the mirror. Since we want to compare the results we have found in the previous section for 
rotating AdS black holes and instability of rotating black holes surrounded by a reflecting mirror found before in the 
four dimensional case [19], and for generic dimensions in [26], here we investigate the stability of Myers-Perry black 
holes under black hole bomb mechanism. To use the asymptotic matching technique, we will need the near and far 
region solutions of the Klein-Gordon equation. 


C. Near-region and far region solutions 


Since near the horizon the Myers-Perry and rotating AdS black holes, in the slow rotation low frequency limit, have 
similar properties, the Klein-Gordon equation near the horizon has the same form (18). Thus, following similar steps, 
the near horizon solution and its far region extension also have the same forms given in (31) and (35). The difference 
will be at the far region solution in which we will now solve in detail. 

In the far region, r — r+ M, since the effects of the black hole can be neglected we have a ~ 0, M ~ 0, A ~ r 2 . 
Therefore, the radial equation is reduced to 


d 2 R{r) i n + 2 dR(r) 
dr 2 


+ 


dr 


2 Ajim 
id — —^ 7 ;— 


R(r) = 0. 


(81) 


The general solution of this differential equation can be expressed in the form of a linear combination of Bessel 
functions of the first and the second kind J and Y as 


R(r) = 


1 


1 + 72 

T 2 


a J ;+ n±i(wr) + BY, . n+i (cor) 


or equivalently in terms of Hankel functions 


R{r) = 


1 


DlH \+i+2L ( Wr ) + D * H l+l±lt ( Wr ) 


(82) 


(83) 


by considering the identities relating Bessel functions to Hankel functions [62]. Here, the relations between the 
integration constants of both forms of the solutions are given by 


Di 



a + ifd 
2 


(84) 


In order to match this far region solution with the near region solution at intermediate regions, we will need the 
small r behavior of the far region solution, which can be easily found from (82) by using the asymptotic forms of the 
Bessel functions at small values [62]. Then, the small r limit of the far region solution becomes 


R{r) 




(1+n+Z) 



a 


P(/ + i±^ + !) 


(85) 


We will also need r —> oo behavior of the far region solution, which can be easily derived by considering (83) form 
of the solution with the asymptotic behavior of Hankel functions for large values of their arguments as given in [62]. 
Then the asymptotic form of the far region solution becomes 


R(r 



Die 




(l + r»)j 


-f) 


D2e 




(!+»)» 


-f) 


( 86 ) 


D. Matching near and far region solutions 


The matching of the near region and far region solutions requires the powers of r of the distinct parts to solutions 
to be the same for both solutions. When MCr-qCi, the near region solution (35) as r —> oo and the far region 
solution (85) as r —> 0 overlap and this matching yields 

r(T^ + i)r(rk + i-2id;) 

r ( ; + 7T + x ) r + l) T (1 - 2iCo)' 


A 


= (!) 


i+± 


(r 


1+n 


_ r l+«)T+^ 


a 


(87) 



14 


and 


P_ 

a 


— 7 r 


/ W \ 2(+n+l 

V2/ 


/ i+n i+n\I+-rr- 

(r + — ) !+« 


r(/ + i±^ + i)r(-i±|) 


r(i + i)r(-^-i)r(i + i-2i) 
r (ifs +1) r ( l + ^ ?(-&-**) 


( 88 ) 


E. Scalar field absorption and decay rates 

Let us calculate here the absorption cross sections and decay rates of a massless scalar field from a rotating black 
hole with single spin. Note that these were investigated before especially in the context of brane world picture in [67]. 


1. The flux and the absorption cross section of the scalar field 
We can calculate the flux using the formula 


J = 


2 i 


dR R dR* 


(89) 


where * denotes complex conjugation. The incoming and outgoing fluxes at radial infinity are calculated as 

Jin =- 1 J^2 1 2 ) 

n 


Jout — \dJl I 

71 


(90) 


The absorbed flux from the black hole horizon located at z = 0 can be calculated from the near region solution of the 
form (31) using the appropriate formula 


J = 


(n + l)(a;+ — x_) z 
2 i 


dR dR* 


which yields 

Jabs = -(n+ l)(a;+ - X-)Cj\A\ 2 = -(w - mfi/ l ) r + +2 |+ 2 , 
where we have used the definition of u> given in (24) at the last step. The absorption probability is given by 

,2 Jabs _r> \ 7rr + +2 


1 — IS 1 ;! 2 = -y— = (u — mQ. h ) 


12 - 


(91) 

(92) 

(93) 


Jin ' ' 2 \D2\ 

Now we are ready to calculate the absorption cross section for a massless scalar field from afl=4 + n dimensional 
singly rotating black hole using the formula [ 68 ] 


2 1 +" 71^ 


m = 


CO 


2+n 


1 + n 


1 + 


1 + n\ (l + 

l 


Using the expression (93) above, we have 


2 n 7r^r 

C 7 ; = (w - mQ h ) + + —;rv— T 


UJ 


2+n 


1 + n 


l + 


(1 - |5z| 2 ) . 

1 + n\ (l + n)\ \A\ 2 


nW. \D 2 1 2 ' 


(94) 


(95) 


Since for small frequency limit (3 <C a, we can approximate |^4| 2 /|Z? 2 1 2 ~ 4|+a| 2 , and using (87), the absorption 
cross section becomes 

+31, ,21-1 2+n, 1+n 1 +„^ « f, . 1 + n\ (l + re)! 


<ti = (w- mCt h ) 


x T 


1 + re 


l (r L + +n - r L _ +n )^ ^ 

(ihl + l) r (iTH + 1 ~ 2iCb ) 


n\l\ 


r(/ + i±2 + i)r( T ^ + i)r(i-2^) 


(96) 


Thus, we see that in the case of superradiance, ui — < 0, the absorption cross section becomes negative implying 

the amplification of the scalar field or equivalently the energy extraction from the black hole. 
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2. Decay rates 


Now let us calculate the decay rate of this scalar field using the formula 


Ti 


_ m _ 

e /3 H (u—mZl h ) _ X ’ 


(97) 


where /3h is the inverse Hawking temperature of the black hole and in the slow rotation limit it is given by (remember 
that r + = ry and r_ = 0 for D > 4) 


Ph 


4tt r 2 + n 
1 + n r++" - ri+" ‘ 


(98) 


The general equation of this expression can be obtained by replacing <7/ in (97) with its corresponding expression (96). 
Now let us calculate its special values when the ratio is integer or half integer. 

When the ratio = k is an integer, the decay rate becomes 


r, = 


9 /_i 3 +n 

ZL 1 7T 2 


(ui — milh) w 
e /3 H (uj—mQh) _ X 2 2l ~ 1 ,+ 


„2+n/„l+ri 


(r L + n - r_ n ) 


k\ 

(2fc)!j 


x [l T 


l + n\(l + n)\ T(i±2) 


nW. r 2 (l + 


+1) M 


no' 2 + 4 - 2 )- 


(99) 


This corresponds to the decay rate of a bosonic field. For n = 0 this expression reduces to the decay rate of the Ken- 
black hole. It is also compatible with the decay rates of five dimensional black holes if the orbital quantum number 
of the scalar field is even [25, 69, 70]. These bosonic terms arise in all dimensions whenever 1/(1 + n) is an integer. 

When the ratio j ^ = k + 1/2 is half integer where k is again an integer, which is only possible for D > 4 and for 
odd D , we can also expand the gamma functions in the general expression of decay rate. Following the similar steps, 
we find that 


r« = 


UJ 


21-1 


5 +n 

7T 2 


e P H (ui-mQ h ) _|_ X 2 2 0 + fc + 1 ) + 

F t i+Iil T / i 

_ 1 ix L _TT (j-l 

T 2 (l + ^ + l) AA [V 2 


(r L + n - r_ n ) T +" +i (n + 1 )[l + 


1 + n\ (l + n)\ 
nlll 


(2k + 1)!! 
(2k + 1)! 


12 


■ 4 U! 


( 100 ) 


This decay rate is similar to the decay rate of a fermionic field. These expressions agree for the decay rates for 
five dimensional black holes [25, 69, 70]. In five dimensions since we can have only integer or half integer values of 
1/(1 + n) the decay rates are either bosonic for even l or fermionic for odd l. For dimensions greater than 5, however, 
1/(1 + n) is not an integer or half integer except for some special cases, the decay rates are not either purely fermionic 
or purely bosonic. These bosonic and fermionic terms, which imply the existence of a conformal symmetry, are also 
obtained in the discussion about high D limit of general relativity in [71]. The possibility of whether the general 
decay rate expression (97) corresponds to a combination of both factors requires further investigations. 


F. Black hole bomb and superradiance instability 

Now we surround the D dimensional singly rotating black hole with an artificial mirror having reflective walls as 
done in [19, 26]. The radius of this wall must be large in order to use the approximations we have adopted. At the 
location of the mirror, r = ro, the radial part of the wave function must vanish, i.e. 


R(r 0 ) = 0. 


( 101 ) 


Considering this condition on the far solution (82), we have 

P _ Jl+(n+ l)/ 2 ( w? ’o) 
a Y l+(n+ i )/2 (ur Q y 


(102) 


where the left hand side of this equation is given in (88). 
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When there is no black hole, the scalar field will develop stable modes which may be labeled by u> = The 
interaction of the scalar field with the black hole will affect this configuration and the frequency of the scalar field 
will be changed. We assume that, in the limits we consider, this change will be small. Thus, we consider that in the 
presence of the black hole the correction to the frequency of these modes will be of the form 

u> = wo + iS, (103) 

where the absolute value of the correction term 6 is assumed to be much smaller than wo- Our aim is to determine 
the correction term iS. Under the small frequency limit, since left-hand side of the ( 88 ) is proportional to co 2l+n+1 , 
we can take in the first approximation that the ratio f3/a in ( 88 ) is vanishing. By considering (102), for a given to 
and ro, we have 


J ;+ n+i(wor 0 ) = 0, (104) 

which has solutions 

Worn = ji + n±i }N , (105) 

where . n+i N are the roots of the equation (104). Their values can be found in [62] or can be obtained by using an 
analytical computer program. In this approximation, the frequency of the scalar field becomes 


uj = wq + iS 


i;+ 2 ±i,jv + 

ro 


(106) 


Using the Taylor expansion in these approximations we have J l+ ^+±(uiro) — iSJ' l+ n±i(ji + 2 ±i jy) in (102) where 
S = 6/r o- Considering all these, from (102) we find that 


- . ( w o r o) ^ 1 i+n ) T (1 + — 2iCj) 

W 7 j; +!f Kro) r(-i±i) r(- T ^^2iw) ’ 

where the factor 

1 -I- 

= (<M0\ 2l+n+1 rq + iir) (r 1 + +n -rl +n ) + ^ 

7 ' 2 / r(i + T ^ : )r(i +1 + i±«)r(z + i±s) ’ 


(107) 


(108) 


is positive. Using the expansion of gamma functions (50) first given by us in [27] and the expansion(52) given in [47], 
5 becomes 


where 


8 = -K 


sinh — i tan 


(l + e)i 


1 + n 


cosh (2ttlo) 


r 7 

/C= 2 ) 


( 1 + ^) 

(2+T^) 


Y i+ n+i(uJoro) 


J' [{ „ + i (cuorp) 


r(i 


1 + n 


— 2ioj) 


(109) 


(HO) 


Thus the perturbation term <5, hence the frequency of scalar field, has both real and imaginary parts. In the general 
case, i.e., if none of these terms vanish, if the frequency of the scalar field satisfies the superradiance condition, and 
if the imaginary part is positive, both an instability and frequency shifts occur for a black hole surrounded by a 
reflective mirror. This is in accordance with the previous works for generic dimensions [26, 27] and also for rotating 
charged black holes in five dimensions [25, 28]. Since the frequency is given by u> = u>o + iS, the real and imaginary 
parts of the frequency are 


Re[w] 


KL 

wo ~ Im[(5] = wq -tan 

ro 


(l + e)n 
1+n 


cosh (27 tw) , 


Im[w] = Re[<5] 


-sinh (27 tw). 

ro 


( 111 ) 


( 112 ) 
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The vanishing of the imaginary part (112) implies that the scalar field does not show an instability and the perturbation 
is oscillating with time. As it is obvious, this is not the case and instability is always present if the superradiance 
condition is satisfied. Note that one can expand the T functions in (109) for some special cases of 1/(1 + n) as in the 
AdS case, without altering the conclusions holding in the general case. 

The time scale of the instability is r = 1/Im[w] = 1/Re[<5]. Since one can deduce the relation <5 ~ l/r^' l+1 ^ +n from 
the above expressions we see that the time scale increases with increasing mirror radius and also increasing spacetime 
dimensions, since we have 


7 - ~ r 2{l+1)+n 
T '0 


(113) 


This implies that in the large D limit, the instability becomes ineffective. 

The real part of the frequency of the waves is inversely proportional to the mirror radius w ~ 1/ro, implying a 
minimum mirror radius to instability to take place. The end point of the instability is when the system reaches the 
critical frequency condition u> c = resulting in a black hole with smaller energy and angular momentum. 


[1] R. Brito, V. Cardoso, and P. Pani, arXiv:1501.06570vl. 

[2] Ya. B. Zel’dovich, Pis’ma Zli. Eksp. Teor. Fiz. 14, 270 (1971), [JETP Lett. 14, 180 (1971)]. 

[3] J. D. Bekenstein, Phys. Rev. D 7, 949 (1973). 

[4] J. M. Bardeen, W. H. Press, and S. A. Teukolsky, Astrophys. J. 178, 347 (1972). 

[5] A. A. Starobinsky, Zh. Eksp. Teor. Fiz. 64, 48 (1973), [Sov. Phys. JETP 37, 28 (1973)]. 

[6] A. A. Starobinsky and S. M. Churilov, Zh. Eksp. Teor. Fiz.65, 3 (1973),[Sov. Phys. JETP 38, 1 (1973)]. 

[7] J. G. Rosa, arXiv:1501.07605vl. 

[8] T. Regge and J. A. Wheeler, Phys. Rev. 108, 1063 (1957). 

[9] C. V. Vishveshwara, Phys. Rev. D 1, 2870 (1970). 

[10] W. H. Press and S. A. Teukolsky, Astrophys. J. 185, 649 (1973). 

[11] W. H. Press and S. A. Teukolsky, Astrophys. J. 193, 443 (1974). 

[12] T. Damour, N. Deruellc, and R. Ruffini, Lett. Nuovo Cimento 15, 257 (1976). 

[13] T. J. M. Zouros and D. M. Eardley, Ann. Phys. (N.Y.) 118, 139 (1979). 

[14] S. L. Detweiler, Phys. Rev. D 22, 2323 (1980). 

[15] H. Furuhashi and Y. Nambu, Prog. Theor. Phys. 112, 983 (2004). 

[16] S. R. Dolan, Phys. Rev. D 76, 084001 (2007). 

[17] Ya. B. Zel’dovich, Zh. Eksp. Teor. Fiz 62, 2076 (1972),[Sov. Phys. JETP 35, 1085 (1972)]. 

[18] W. H. Press and S. A. Teukolsky, Nature (London) 238, 211 (1972). 

[19] V. Cardoso, O. J. C. Dias, J. P. S. Lemos, and S. Yoshida, Phys. Rev. D 70, 044039 (2004). 

[20] M. J. Strafuss and G. Khanna, Phys. Rev. D 71, 024034 (2005). 

[21] S. Hod and O. Hod, Phys. Rev. D 81, 061502 (2010). 

[22] J. G. Rosa, J. High Energy Phys. 06 (2010) 015. 

[23] H. Witek, V. Cardoso, C. Herdeiro, A. Nerozzi, U. Sperhake, and M. Zilhao, Phys. Rev. D 82, 104037 (2010). 

[24] S. R. Dolan, Phys. Rev. D 87, 124026 (2013). 

[25] A. N. Aliev, and O. Delice, Superradiance, scalar field absorption and black hole bomb for five dimensional charged rotating 
black holes, (unpublished) (2009). 

[26] J. P. Lee, J. High Energy Phys. 01 (2012) 091. 

[27] T. Durgut, Ms. thesis, Marmara University Physics Department, 2013. 

[28] A. N. Aliev, J. Cosmol. Astropart. Phys. 1411(2014)029 

[29] J. C. Degollado, C. A. R. Herdeiro, and H. F. Runarsson, Phys. Rev. D 88, 063003 (2013). 

[30] J. C. Degollado and C. A. R. Herdeiro, Phys. Rev. D 89, 063005 (2014). 

[31] S. Hod, Phys. Rev. D 88, 064055 (2013). 

[32] R. Li, J. Zhao, and Y. Zhang, Conimun. Theor. Phys. 63, 569 (2015). 

[33] R. Li and J. Zhao, Phys. Lett B 740, 317 (2015). 

[34] J. M. Maldacena, Adv. Theor. Math. Phys. 2, 231 (1998). 

[35] S. W. Hawking and H. S. Reall, Phys. Rev. D. 61, 024014 (1999). 

[36] V. Cardoso and O. J. C. Dias, Phys. Rev. D 70, 084011 (2004). 

[37] N. Uchikata, S. Yoshida, and T. Futamase, Phys. Rev. D 80, 084020 (2009). 

[38] S. S. Gubser, Phys. Rev. D 78, 065034, (2008). 

[39] A. N. Aliev and O. Delice, Phys. Rev. D 79, 024013 (2009). 

[40] O. J. C. Dias, P. Figueras, S. Minwalla, P. Mitra, R. Monteiro, and J. E. Santos, J. High Energy Phys. 08 (2012) 117. 

[41] V. Cardoso, O. J. C. Dias, and S. Yoshida, Phys. Rev. D 74, 044008 (2006). 

[42] H. K. Kunduri, J. Lucietti, and H. S. Reall, Phys. Rev. D 74, 084021 (2006). 

[43] H. Kodama, Prog. Theor. Phys. Suppl. 172, 11 (2008). 

[44] K. Murata, Prog. Theor. Phys. 121, 1099 (2009). 



18 


[45] H. Kodama, R. A. Konoplya, and A. Zliidenko, Phys. Rev. D 79, 044003 (2009). 

[46] V. Cardoso, O. J. C. Dias, G. S. Hartnett, L. Lehner, and J. E. Santos, J. High Energy Phys. 04 (2014) 183. 

[47] M. Wang and C. Herdeiro, Phys. Rev. D 89, 084062 (2014). 

[48] B. Carter, Commun. Math. Phys, 10, 280 (1968). 

[49] G. W. Gibbons, H. Lu, D. N. Page, and C. N. Pope, J. Geom. Phys. 53, 49 (2005). 

[50] G. W. Gibbons, H. Lu, D. N. Page, and C. N. Pope, Phys. Rev. Lett. 93, 171102 (2004). 

[51] S. W. Hawking, C. J. Hunter, and M. M. Taylor-Robinson, Phys. Rev. D 59, 064005 (1999). 

[52] G. W. Gibbons, M. J. Perry, and C. N. Pope, Class. Quantum Grav. 22, 1503 (2005). 

[53] M. M. Caldarelli, G. Cognola, and D. Klemm, Class. Quantum Grav. 17, 399 (2000). 

[54] P. T. Chrusciel, D. Maerten, and P. Tod, J. High Energy Phys. 11 (2006) 084. 

[55] O. J. C. Dias, P. Figueras, R. Monteiro, and J. E. Santos, J. High Energy Phys. 12 (2010) 067. 

[56] B. Gwak and B. H. Lee, Phys. Rev. D 91, 064020 (2015). 

[57] H. K. Kunduri and J. Lucietti, Phys. Rev. D 71, 104021 (2005). 

[58] M. Vasudevan and K. A. Stevens, Phys. Rev. D 72, 124008 (2005). 

[59] C. Muller, Spherical Harmonics, Lecture Notes in Mathematics, (Springer, New York, 1966). 

[60] H. Suzuki, E. Takasugi, and H. Umetsu, Prog. Theor. Phys, 100, 491 (1998). 

[61] H. T. Cho, A. S. Cornell, J. Doukas, and W. Naylor, Phys. Rev. D 80, 064022, (2009). 

[62] M. Abramowitz and A. Stegun , Handbook of Mathematical Functions (Dover Publications, New York, 1970). 

[63] R. Emparan, and R. C. Myers, J. High Energy Phys. 09 (2003) 25. 

[64] O. J. C. Dias, P. Figueras, R. Monteiro, J. E. Santos, and R. Emparan, Phys. Rev. D 80, 111701 (2009). 

[65] R. Emparan and H. S. Reall, Living Rev. Rel. 11, 6 (2008). 

[66] E. Berti, V. Cardoso, and M. Casals, Phys. Rev. D 73, 024013 (2006). 

[67] S. Creek, O. Efthimiou, P. Kanti, and K. Tamvakis, Phys. Lett. B 656, 102 (2007). 

[68] S. S. Gubser, Phys. Rev. D 56, 4984 (1997). 

[69] J. M. Maldacena and A. Strominger, Phys. Rev. D 56, 4975 (1997). 

[70] O. J. C. Dias, R. Emparan, and A. Maccarrone, Phys. Rev. D 77, 064018 (2008). 

[71] R. Emparan, S. Suzuki, and K. Tanabe, J. High Energy Phys. 06 (2013) 009. 



